Correlation Effects on Magnetic Anisotropy in Fe and Ni 
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We calculate magnetic anisotropy energy of Fe and Ni by taking into account the effects of strong 
electronic correlations, spin-orbit coupling, and non-coUinearity of intra-atomic magnetization. The 
LDA+U method is used and its equivalence to dynamical mean-field theory in the static limit is 
derived. The effects of strong correlations are studied along several paths in {U, J) parameter space. 
Both experimental magnitude of MAE and direction of magnetization are predicted correctly near 
U = 1.9 eV, J = 1.2 eV for Ni and U = 1.2 eV, J = 0.8 eV for Fe. The modified one-electron 
spectra by strong correlations are emphasized in conjunction with magnetic anisotropy. 
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I. INTRODUCTION 



One of the long-standing problems that are still short 
of detailed understanding is to explain the magneto- 
crystalline anisotropy energy of magnetic materials con- 
taining transition-metal elements, especially that of Fe, 
Co, and Ni. Magnetic anisotropy is the dependency of 
internal energy on the direction of spontaneous magneti- 
zation. Generally, the magnetic anisotropy energy term 
possesses the symmetry of the crystal. It is therefore 
called magneto-crystalline anisotropy or crystal magnetic 
anisotropy. In transition metals, most of the magnetic 
moment conies from spin polarization. Within the non- 
relativistic approach, there is no term coupling to the 
spin degrees of freedom. The magnetic moment, there- 
fore, can points to an arbitrary direction. The mag- 
netic anisotropy energy is an relativistic phenomenon of 
spin-orbit coupling. With spin-orbit coupling, the spin 
degrees of freedom interact with the spatial anisotropy 
through the coupling to the orbital degrees of freedom. 
This induces a preferred direction of spins. This hap- 
pens even though the total angular momentum is com- 
pletely quenched. Since spin-orbit coupling couples indi- 
vidual spin degrees of freedom to individual orbital de- 
grees of freedom, the fact that total orbital magnetic mo- 
ment is quenched does not affect the visibility of spatial 
anisotropy to magnetic moment. 

The primary difficulty toward investigating MAE has 
been attributed to the fact that the MAE of the met- 
als are as small as of the order of 1 /ieV^/atom. At 
low temperature (T — A.2K), MAE is of the order of 
60 ^eV/aiom for hep Co, 2.8/iey/atom for fee Ni, and 
1.4 /xe^/atom for bcc Fc. With advances in the accurate 
total energy method combined with the development of 
faster computers, attempts have been made to calculate 
MAE from first-principles for bulk crystalline Fe, Co, 
and Ni. While the magnitude of the MAE have been 
predicted for all the three metals, the correct easy axis 



for Ni has not been predicted so far. Here we present 
new results that predict the correct easy axis for Ni as 
well as reproduce the previous results in relevant lim- 
its. This result is obtained by incorporating spin-orbit 
coupling, non-coUinear nature of of intra-atomic magne- 
tization, and strong correlation effects. We also suggest 
why the previous approaches have failed in obtaining the 
correct easy axis for Ni. 

Experimentally, there are various ways of measuring 
the magnetic anisotropy energy. These range from torque 
measurements, saturation magnetic fields and ferromag- 
netic resonance methods. However, torque measurements 
are more often used and are considered to be reason- 
ably precise. The main principle behind the torque ex- 
periments is that in the presence of an external mag- 
netic field misaligned with respect to the easy axis (or 
equivalently, direction of the crystal magnetization), the 
magnetic moment experiences a torque. In equilibrium, 
this torque is balanced by the crystal anisotropy torque. 
Torque magneto-meters are used to precisely measure 
this torque. Since the crystal anisotropy torque is just 
the angular derivative of the magnetic anisotropy energy, 
a measurement of the torque gives an indirect estimate of 
the anisotropy energy. The torque measured at various 
angles are interpolated to a corresponding analytic ex- 
pression such as angular derivatives of total energy. The 
analytic expression of total energy can be expressed in 
terms of parameters based on crystal symmetry. 

The simplest form of crystal magnetic anisotropy is 
uniaxial anisotropy, for example, in hexagonal cobalt 
with easy direction parallel to the c axis of the crystal at 
room temperature. As the internal magnetization rotates 
away from the c axis, the anisotropy energy increases 
with increase of (/>, where 4> is the angle between the c 
axis and the internal magnetization. We can expand this 
energy in a series of powers of sin^ 0: 



Ea = Ki sin' 



K2 sin'' 



(1) 



where Ki and K2 are constants. 



2 



For cubic crystals such as iron and nickel, the 
anisotropy energy can be expressed in terms of the di- 
rectional cosines (01,02,013) of the internal magnetiza- 
tion with respect to the three cubic edges. Because of 
the high symmetry of the cubic crystal, the anisotropy 
energy can be expressed in a fairly simple way: 

Ea = Ki {ajaj + ajal + alaf) + K2 (01O2O3) H , 

(2) 

where Ki and K2 are constants. 

The constants Kn obtained by interpolating the torque 
are substituted back to get the expression for the total 
energy. Magnetic anisotropy energy is calculated by sub- 
tracting energies at two different direction of magnetic 
moments, e.g. [001] and [111] 

To measure the torque, specimens are formed into 
short cylinders. The specimens are mounted in a car- 
riage held by two torsion fibers that are fastened to a 
rigid support at the top and to a circular scale S at the 
bottom. When the field is excited in the electromagnet, 
the crystal tends to turn so that the direction of easy 
magnetization is parallel to the field. The torque so pro- 
duced is balanced by turning the bottom of the lower 
fiber until the crystal regains its original orientation as 
determined by reflection of a light beam from the mirror. 
The scale reading S2 is then compared with the original 
reading Si with H — 0. 6*2 — 6*1 is a measure of the 
torque. The orientation of the crystal axes with respect 
to the applied field is varied by turning the electromag- 
net, which is mounted on a heavy bearing, and noting 
its position on a suitable scale, S'. One then plots the 
torque against the crystal orientation and deduces from 
the curve the crystal anisotropy constant Kn- 

Early attempts to explain magnetic anisotropy are 
based on the interaction between the magnetization and 
the lattice through spin-ort»it coupling combined with 
band theoryld. Later, Brooka3 used an itinerant electron 
model and orbital angular momentum quenching in cubic 
crystals to explain MAE. Treating the spin-orbit coupling 
as a perturbation, the nontrivial magnetic anisotropy 
came out at the fourth-order. The correct directions of 
magnetization were obtained for Fe and Ni. The differ- 
ence between easy axes in these metals was attributed to 
the different lattice structures of them. Based on Ni-Fe 
alloy data, however, a much closer correlation between 
the anisotropy and the number of valence electrons was 
observed. What is the factor that determines the easy 
axis has been anothst pVinanswered question since then. 
In subsequent paperg3'QH the calculations became more 
and more refined. Finally the discussion centered around 
the importance of degenerate states along symmetry lines 
in the Brillouine zone and plausible explanalipns of the 
origin of magnetic anisotropy were providedElQ. 

Contemporary studies of magnetic anisotropy energy 
have centered on firstjprinciples calculations, using den- 
sity functional theorylj within tlpje framework of local spin 
density-approximation (LSDA)tl. Eckardt, Fritschc, and 
Noffkety were able to get the values of the right order of 
magnitude for Fe and Ni, but with incorrect easy axis for 



Fe. Daaiderop, Kelly, and SchuurmansHl used the force 
theoremtd to obtain the correct order of magnitude of the 
MAE for Fe, Co, and Ni, but incorrect easy axes were 
found for Co, and Ni. They also observed that chang- 
ing the number of valence electrons, hence removing the 
LDA artificial X2 pocket would restore the correct easy 
axis. This observation is similar to Brooks' observation 
of the close correlation bewteen the direction of magnetic 
moment and the number of valence electrons. Wang, Wu, 
and Freeman developed state tracking method to study 
magnetic Mjsotropy of Fal3. Trygg, Johansson, Eriksson, 
and WillisEj improved the method while using fully self- 
consistent approach. Orbital polarization|-was also in- 
corporated, which was suggested by Jansenlla. While the 
correct easy axis was obtained for Co and Fe, in case of Ni 
the calculation still gave the wrong easy axis. This is the 
best result before the current work. Schneider, Erickson, 
and Jansent£l used torque instead of energy difference 
to obtain the same result as that of Trygg, Johansson, 
Eriksson, and Willis. They also treated the spin-orbit 
coupling constant as an adjustable parameter. They suc- 
ceeded in restoring the correct easy axis for Ni with an 
unphysically large value of spin-orbit coupling. This work 
suggested a close relation between magnetic anisotropy 
energy and the strength of spin-orbit coupling. Halilov, 
Perlov, Oppeneer, Yaresko ,and Antonov also scaled spin- 
orbit coupling in order to enlaxge its effect on the MAE 
and confirmed previous resultal3. 

We believe that the physics of transition metal com- 
pounds is intermediate between atomic limit where the 
localized d electrons are treated in the real space and fully 
itinerant limit when the electrons are described by band 
theory in k space. A many-body method incorporating 
these two impottant limits is the dynamical mean-field 
theory (DMFTp. The DMFT approach has been ex- 
tensively used to study model Hamiltonian of correlated 
electron systems in the weak, strong and intermediate 
coupling regimes. It has been very successful in describ- 
ing the physics of realistic systems, like the transition 
metal oxides and, therefore, is expected to treat properly 
the materials with d or / electrons. 

We take a new view that the correlation effects within 
the d shell are important for the magnetic anisotropy 
of 3d transition metals like Ni. These effects are not 
captured by the LDA but are described by Hubbard-like 
interactions presented in these systems and need to be 
treated by an extension of first principles mattods such 
as DMFT. Since, DMFT reduces to LDA-hUO in static 
limit, we adopt LDA-I-U method to attack the problem 
of magnetic anisotropy of 3d transition metals. 

Another effect which has not been investigated in the 
context of magnetic anisotropy calculations is ±lie non- 
collinear nature of intra-atomic magnctizationtj. It is 
expected to be important when spin-orbit coupling and 
correlation effects come into play together. We show that 
when we include these new ingredients into the calcula- 
tion we solve the long-standing problem of predicting the 
correct easy axis of Ni. Part of this work has been pub- 
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lished elsewhere0. 

The remainder of this paper is organized as follows. 
Section |l| presents DMFT extended to relativistic elec- 
tron structure problem, its equivalence to LDA+U in 
static limit. Section |l| describes calculated MAE. Sec- 
tion [V discusses the calculated MAE in conjunction with 
band structure. Section ^ presents how to include other 
relativistic corrections than spin-orbit coupling and dis- 
cuss the effects of spin-other-orbit coupling. Section VI 
is a summary and conclusion. 



II. DMFT AND LDA+U 

The LDA method is very successful in many materials 
for which the one-electron model of solids works. How- 
ever, in correlated electron system this is not always the 
case. In strongly correlated situations, the total energy 
is not very sensitive to the potential since the electrons 
are localized due to the interactions themselves, and the 
lack of sensitivity of the functional to the density, does 
not permit to devise good approximations to the exact 
functional in this regime. For example, when the Mott 
transition takes place the invertibility condition is not 
satisfied. Our view is that these difficulties cannot be 
remedied by using more complicated exchange and cor- 
relation functionals in density functional theory. DMFT 
is a me|tiiod successfully describing strongly correlated 
systemstiand has been extended to electronic structure 
problem£3. In this work we utilize the new functional 
formulation of the electronic structure problem and ex- 
tend it to relativistic case. 

The basic idea is to introduce another relevant variable 
in addition to the density p and the magnetic moment 
density m, namely the local Green function. The lat- 
ter is defined by projecting the full Green function onto 
a separate subset of correlated "heavy" orbitals distin- 
guished by the orbital index a and the spin index a from 
a complete set of orbitals xZi^^ — R) = XaR ^ tight- 
binding representation which we assume for simplicity 
to be orthogonal. The local Green-imction is therefore 
given by a matrix G with elementsES 



Gl^ {lu:, R) = - {^clj,{toj)cl+{ioj)) = (3) 
- j xr(r-R) (^(r,zc.)^+(r',zc.))x^'(r'-R)drdr'. 

We then construct a functional r[p, 111,(5] which gives 
the exact free energy at a stationary point. 

To describe the new method, we consider fermion 
system under an external potential T4xt and an ex- 
ternal magnetic field h. For relativistic effects spin- 
orbit coupling, whose effects are important for magnetic 
anisotropy calculations, are also considered. Spin-orbit 
couplinffi is included according to the suggestion by An- 
dersenCJ. It is useful to introduce the notion of the Kohn- 
Sham potential VkSj the Kohn-Sham magnetic field Hks 
and its dynamical analog («Wn)- They are defined 



as the functions that one needs to add to the kinetic en- 
ergy matrix so as to obtain a given density and spectral 
function of the heavy orbitals namely: 

p(r)= r^Tr,(rs|[(zc^„ + VV2- Vks)I 



- 2AiBS-hKs-e(r)l-s-S]]-i|rs)e 
m(r) = -2/isT^Tr,(rs|s[(zc^„ + VV2-l^Ks)I 

-2fiBS ■ hKs - e(r)l • s - I rs) e''^"°^(5) 

where Tig is the trace over spin space, 1 and s are one- 
electron orbital and spin angular momentum operator, 
respectively. The spin angular momentum operator is 
expressed in terms of Pauli matrices s = a/2 and I is 
2x2 unit matrix. Vks and Hks are functions of r. The 
chemical potential fj, is set to zero throughout the current 
section, and S is given by 

EEES(r,r',zc.)= J2 Xa* (r ^ '^)Kb' MXh' (r' -R.) . 



ahaa' R 



(6) 

^(r) determines the strength of spin-orbit coupling and 
in practice is determinedB^ by radial derivative of the 
I = component of the Kohn-Sham potential inside an 
atomic sphere: 



2 dV^sjr) 
c? dr 



(7) 



When spin-orbit coupling is present, the intra-atomic 
magnetization m(r) is not coUinear, i.e., the direction of 
magnetization depends on the position r. Therefore, the 
magnetization must be treated as a general vector field, 
which realizes non-collinear intra-atomic nature of this 
quantity. Such getieral magnetization scheme has been 
recently discussedEy 

In terms of these quantities and the matrix of local 
interactions U, we write down the DMFT+LSDA func- 
tional: 

Tlsda+dmptIp, Vks, hKs, G, S) = 

J2 e"°^ Tr log[(zc^ + - Vks)1 - 2/iss • hxs 

-^(r)l • s - S] - y dvVKs{r)p{r) + J drm(r) • hKs(r) 

-^e'"°^Tr[SMG(zc^)] + ( dvVe^t{r)p{r) 

- /drh(r).m(r) + i 

J 2 J \r-r'\ 

+E^rHp, m] + J2mG] - <i>nc]- (8) 

R 

^[G] is the sum of the two-particle irreducible local di- 
agrams constructed with the local interaction matrix U, 
and the local heavy propagator G. $dc is the so-called 



4 



double counting term which subtracts the average energy 
of the heavy level already described by LDA. Expression 
ensures that the Greens function obtained from its 
extremization will satisfy the Luttinger theorem. 

E^^^^Ip, m] is the LSDA exchange correlation energy. 
Since the exact exchange correlation energy functional is 
not known, the usefulness of this approach is due to the 
existence of successful approximations to the exchange 
correlation energy functional as Kohn and Sham pro- 
posed. When nontrivial magnetic moment is present, the 
exchange correlation energy functional is assumed to be 
dependent on density and magnetization: 



(9) 



dr/xc[/3(r),TO(r)]m(r)(r), (10) 



where m — |m|. 

The functional (||) can be viewed as a functional of 
six independent variables, since the stationary condition 
in the conjugate fields reproduces the definition of the 
dynamical potential and the Weiss field. Extremizing it 
with respect Vks j hxs i and S lead us to compute the den- 
sity (Eq. (P), the magnetic moment density (Eq. (||)), and 
the Green function G^^ {iui) (Eq. (p^)), respectively. The 
Kohn-Sham potential Vks (r) and Kohn-Sham magnetic 
field hKs(r) are obtained by extremizing the functional 
with respect to p(r) and m(r): 



VKs(r) = Vext(r)+ Jdr' 



, p(r') , 5E^^^[pM 



5p{v) 



,11) 



hKs(r) = h(r) 



5E^^^[pM 
(5m(r) 



(12) 



Extremizing with respect to G^J^ yields the equation for 
self energy 



(5$ 



(5$ 



DC 



5Gll'{iu) 



(13) 



The physical meaning of the dynamical potential E is 
parallel to the meaning of the original Kohn-Sham po- 
tential Vks- it is the function that one needs to add to 
the correlated block of the one-electron Hamiltonian in 
order to obtain the desired local Green function: 



G^,-'(^c.)=^[^c.-7^'< 



(14) 



where H^^'^ = (x^kK-V^ + Vks)1 + 2pBS ■ Hks + e(r)l • 
s|xbk) is the one-electron Hamiltonian in k-space. It 
is the frequency dependence of the dynamical potential 
which allows us to treat Hubbard bands and quasiparticle 
bands on the same footing. 

In general, an explicit form of $[G] is not available. 
DMFT maps the DMFT-I-LSDA function to an Ander- 
son impurity model. Self-consistency equations obtained 



in this way are used to find the self energy (|13|). In 
this paper we confine ourselves to zero temperature and 
make an additional assumption on solving the impurity 
model using the Hartree~Fock approximation. In this 
limit an explicit form of ^[G] is available and DMFT 
self-consistency loop is unnecessary. We first figure out 
the Coulomb interaction by considering a Hartree-Fock 
averaging of the original expression for the Coulomb in- 
teraction given by 



i^^(aa6a'|^|cada')cr<'+c5'<. (15) 

aa' abed 

In this limit, the sum of local graphs reduce to 

abcda 

+ 2 X] ^^o.^cd - Jabcd)'^ab'^ld 
abcda 

-\Y. JabcdK^^n-r- (16) 



ahcda 



Here, the matrices Uabcd = {ac\vc\bd) and Jabcd 
{ac\vc\dh) have the following definitions: 



Uabcd - j xl%r)xT{r')vc{r-r')xl{v)Xd{^')dvdv', 
Jabcd = /xa(r)"xr(r')«c(r-r')x5(r)xnr')drdr', 



where the Coulomb interaction vc{v — y') has to take into 
account the effects of screening by conduction electrons. 
Note that the matrices Uabcd and Jabcd are spin inde- 
pendent since the Coulomb interaction is independent of 
spin. The occupancy matrix n'^^ is a derived quantity 
of the Green function: 



(17) 



Notice that when spin-orbit coupling is taken into ac- 
count, the occupancy matrix becomes non-diagonal with 
respect to spin index even though the interaction matri- 
ces Uabcd and Jabcd are spin independent. 

The self energy EJ^^ now takes the from for spin diag- 
onal elements 



^ah = X! U<^bcdnJ " -f y^([/afccrf - Jabcd)n'^^ 

cd cd 

and for spin off-diagonal elements it is given by 



(18) 



cd ^^ab (*^) 



(19) 



The off-diagonal elements of the self energy only present 
when spin-orbit coupling is included, hence a relativistic 
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effect. To make it more physically transparent we can 
introduce magnetic moments at the given shell by 



aa' 



(20) 



where /x runs over x,y,z for Cartesian coordinates, or 
over,— 1, 0, +1 (z, ±) for spherical coordinates. Relativis- 
tic correction from strong correlations can be written in 
physically transparent form 



'cd 



9 2^ "'ab Jabcdm^^ 



abcdcr 



ibcd 



abed 



and in principle has room for further generalization of 
exchange matrix Jabcd to be anisotropic, i.e depend on 
ijij'- T^^' 



bed' 



Part of the energy added by ^[G] is already included 
in LSDA functional. The double counting term <i>(jc is 
added to subtract this already included part of ^[G]. It 
was proposedc3 that the form for <I>[G] is 



Model 
dc 



= -Un{n~l)'~-J[n^n^ ~l)+n^{n^ -1)], (22) 



where 
U 

J 

and n"^ 



ab 



(23) 



U 



2/(2/ 



1 



ab\-\ba))X2A) 



J2a<a^ and n 



ab 



The subtraction by 



1 is made to take the self-interaction into account. This 
generates the self energy in the form: 



^ab 



cd 



Uabcdn^^ " + Y^iJJabcd ~ Jabcd)n'^d 



The DMFT self consistency equation identifies the 
Green function of the original model and the Green func- 
tion of the mapped impurity model to find the self energy. 
Now that we can express the sum of local graphs $[G] 
in terms of the original Green function, the DMFT loop 
need not to be performed. The problem is now reduced to 
extremizing the functional [Eq. m}] w ith the expression 
for the sum of local -graphs [Eq. pq)], which is exactly 
the LDA-t-U methodti. 

The DMFT functional and its static correspondent 
LDA-I-U functional are defined once a set of projectors 
{xJJ(r)} and a matrix of interactions Uabcd and Jabcd are 
prescribed. When / orbitals are used as the projection 
operators, the matrix is expressed in terms of Slater pa- 
rameters F''. For a = lm,b = lk,c = l'm',d = I'k' 
and representing ^^(r) = ^/„(r)(l, 0)"^, where (r) = 
(j)i{r)i''Yi„i{f), we can express the matrices Uabcd and 
Jabcd in the following manner: 



min(2/,2/') 

{lml'm'\-\lkl'k') = ^ 



/"=0,2,... 



2/" + l 



(30) 



/ -1 \ra" ^l" m" —ni—k y-il" ra" — k' ~va' 

^ \~^) ^Iklm ^I'm'l'k' 



min(2/.2/') 

{lml'm'\l\l'k'lk) = Yl WTl""^'"' (^'^ 



i"=0,2,. 



/ -1 \m" ^l" 7n" —m — k' /^l" m" —k~ra' 

^ K'^) ^I'k'lrn ^I'm'lk 

where the quantities Cf^, are the Gaunt coefficients 
which are the integrals of the products of three spher- 
ical harmonics 



C'ti,= / YL{v)Yl,{v)YL..{v)dv. 



(32) 



The quantities i^^") and F(^) are given by the following 
radial integrals 



- 5abU{fl~]^)+5abJ{n'' 




(25) 




r r' 


V^^"" \ T 

^ab — ^ 2^ ■Jabcdn^d ■ 
cd 




(26) 




r r'" 
/ J.ll"+1 



—^fir)4ir')drdr' 



(33) 



As an example, when only the effect of U is under 
investigation, the U and J matrices are Uabcd — SabScdU, 
Jabcd = SadScbU, U ~ U, and J = 0. This simple U and 
J matrices make it possible to write down corrections to 
LSDA functional and LSDA Kohn-Sham potential: 



'6a 



-Un 



(27) 



When 1 = 1', the quantites F^") and F^-^^ are equal and 
have a name of Slater integrals which for s-electrons are 
reduced to one constant F^^\ for p-electrons there are 
two constants: F^°\ F^^), for d's: F^, F^^\ F^, etc. 
In this case, the expressions for U and J are reduced to 

{m,m"\vc'\m',m"') =J2''f'{m,m',m",m"')F'', (35) 
k 

where < k < 21, and 



Kb = U{-5ab 



(28) 
(29) 



6 



Slater integrals can be linked to Coulomb and Stoner 
parameters U and J obtained from LSDA supercell pro- 
cedures via U = F° and J = {F"^ + F^)/IA. The ratio 
F'^/F'^ is to a good accuracy a constant ^ 0.625 for d 
electrons. For / electrons, the corresponding expression 
isU = F° and J = (286^2 + 195^" + 250i^6)/g435 

To summarize, we have shown the equivalence of 
Hartree-Fock approximation of DMFT and LDA+U 
method. LDA+U is the method proposed to over- 
come the. difficulties of LDA when strong correlations are 
presenicJ. Since the density uniquely defines the Kohn- 
Sham orbitals, and they in turn, determine the occu- 
pancy matrix of the correlated orbitals, once a choice of 
correlated orbitals is made, we still have a functional of 
the density alone. However it is useful to proceed with 
Eq. (^) , and think of the LDA -I- U functional as a func- 
tional of p, Vks, m, hxsj Gab i ^^'^ ^at ' whose min- 
imum gives better approximations to the ground-state 
energy in strongly correlated situations. Allowing the 
functional to depend on the projection of the Kohn-Sham 
energies onto a given orbital, allows the possibility of or- 
bitally ordered states. This is a major advance over 
LDA in situations where this orbital order is present. As 
recognized many years ago, this is a very efficient way of 
gaining energy in correlated situations, and is realized in 
a wide variety of systems. 



III. NUMERICAL CALCULATION OF MAE 

We calculate MAE by taking the difference of two to- 
tal energies with different directions of magnetization 
[MAE=£'(111) - E{001)]. The total energies are ob- 
tained via fully self consistent solutions. Since the total 
energy calculatkm requires high precision, full potential 
LMTO methodcj has been employed. For the k space 
integrationjve follow the analysis given by Trvgg and 
co-workersLj and use the special point methodcj with 
a Gaussian broadening23 of 15 mRy. The validity and 
convergence of this procedure has been tested in their 
work EJ. For convergence of the total energies within de- 
sired accuracy, about 15000 fc-points per Brillouine zone 
are needed. We used 28000 fc-points to reduce possi- 
ble numerical noise, while the convergency is tested up 
to 84000fc-points. Our calculations include non-spherical 
terms of the charge density and potential both_within 
the atomic spheres and in the interstitial regionEfl. All 
the low-lying semi-core states are treated together with 
the valence states in a common Hamiltonian matrix in 
order to avoid unnecessary uncertainties. These calcu- 
lations are spin polarized and assume the existence of 
long-range magnetic order. Spin-orbit coupling is ipir 
plemented according to the suggestions by AndersencJ. 
We also treat magnetization as a general vector field, 
which realizes non-collinear intra-atomic nature of this 
quantity. Such geaeral magnetization scheme has been 
recently discussedE3. 

We first test our method in case of LDA {U = J = 



0). To compare with previous calculations, we turn 
off the non-collinearity of magnetization which makes 
it coUinear with the quantization axis. The calcula- 
tion gives correct orders of magnitude for both fee Ni 
(0.5 iJ.eV) and bcc Fe (0.5 fJ.eV) but with the wrong easy 
axis for Ni, which is the same result as the previous oneEj. 
Turning on the non-collinearity results in a a larger value 
of the absolute value of the MAE (2.9 p,eV) for Ni but 
the easy axis predicted to be (001) which is still wrong. 
The magnitude of the experimental MAE of Ni is 2^ ^eV 
and the easy axis is aligned along (111) directionEll. For 
Fe, the non-collinearity of magnetization changes neither 
MAE (0.5 neV) nor the easy axis (001) from the colhnear 
result. The magnitude of the experimental MAE of Fe is 
1.4 iJ.eV and the easy axis is aligned along (001) direction. 
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FIG. 1: Ni. Contour plot of magnetic moment as a function of 
U and J. The contour is drawn at 0.018 /^b interval, which is 
2.9% of the experimental value of magnetic moment 0.606 fj,B- 




0.5 1 1.5 2 2.5 



U(eV) 

FIG. 2: Fe. Contour plot of magnetic moment as a function 
of U and J. The contour is drawn at 0.063 fiQ interval, which 
is 2.9% of the experimental value of magnetic moment 2.2 fi^. 

We now describe the effect of correlations, which is 
crucial in predicting the correct axis of Ni. We first study 
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the effects of strong correlations on magnetic moments. 
We have scanned the {U, J) parameter space to obtain 
magnetic moment as a function of U and J (see Figs. 
|l| and H). Magnetic moment increases as U increases, 
but decreases as J increases for both Ni and Fe. The 
magnetic moments of Ni and Fe change up to 20% in the 
parameter range. Notice that weU-defined contours exist 
where the magnetic moments are close to experimental 
values for both Ni and Fe. In comparison, the magnetic 
moment of Fe depends more strongly on J than that, of 
Ni. This result is in agreement with an earlier workE3. 

We now discuss our calculated MAE. The load of com- 
puting MAE is very heavy. Rather than calculating the 
quantities in the ([/, J) parameter space, we follow three 
paths: two paths with constant J, and the path with 
experimental magnetic moments. 




FIG. 3: The magneto-crystalline anisotropy energy MAE = 
£(111) - £(001) for Ni (top) and Fe (bottom) as functions 
of U. The experimental MAEs are marked by arrows for Fe 
(1.4 neV) and Ni (-2.8 /leV). The value of J is held to 0. 

We first walk along a path with J = eV. For Ni 
(see Fig. as U increases, the MAE of Ni smoothly 
increases until U reaches 2.5 eV and then smoothly de- 
creases up to the value 3.8 ^eV. Around U = 3.9 eV, 
the MAE decreases abruptly to negative value. Around 
U — 4.0 eV, the experimental order of magnitude and the 
correct easy axis (111) are restored. The change from the 
wrong easy axis to the correct easy axis occurs over the 
range of SU ^ Q.2eV, which is the order of spin-orbit 
couphng constant O.leV). 

For Fe, the MAE decreases on increasing U to negative 
values, where the magnetization takes the wrong axis. 
From U = 2.7 eV, it increases back to the correct direc- 



tion of easy axis (positive MAE). Around t/ = 3.5 eV, it 
restores the correct easy axis and the experimental value 
of MAE is reproduced. 

In this study on effects of intra-atomic repulsion, we 
see that it is possible to predict the correct magnetic 
anisotropy at a nontrivial U. Notice, however, that the 
calculated magnetic moments are about 20% larger than 
those of experiments for both Ni and Fe. As seen in 
the Figs. 1^ and ^, the magnetic moment decreases as J 
increases. This lead us to follow a path with nontrivial 
J. 
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FIG. 4: The magneto-crystalline anisotropy energy MAE — 
£(111) - £(001) for Ni (top) and Fe (bottom) as functions 
of U. The experimental MAEs are marked by arrows for Fe 
(1.4 fJ,eV) and Ni (—2.8 fieV). The values of J are fixed to 
0.9 for both Ni and Fe. 

We pick J = 0.9 eV that seems to be a canonical value 
of J. For Ni, (see Fig. ^), as U increases, MAE increases 
to 15 ^eV till U = 1.6 eV. Then MAE decreases, chang- 
ing the easy axis from (001) to (111) at C/ = 2.6 eV. The 
experimental value of MAE with the correct easy axis is 
predicted at 2.7 eV. 

For Fe, MAE increases to 10 ^eV tiU U = 1.9 eV. Then 
MAE decreases, changing the easy axis to (111) direction 
at [/ = 2.2 eV. The experimental value of MAE with the 
correct easy axis is predicted at 1.7 eV. 

In this study of effects of intra-atomic repulsion ([/) 
with nontrivial J, we are also able to predict the correct 
magnetic anisotropy at a nontrivial U, much less than 
that of the case with trivial J. Notice that the calculated 
magnetic moments are about 5% larger than those of 
experiments, a big improvement over the case with J = 0. 
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To treat properly the correlation effects on calculated 
anisotropy energy, the intra-atomic repulsion U and ex- 
change J should be taken into account spontaneously. 
We, therefore, follow a path in ([/, J) space determined 
by fixed values of magnetic moments close to experimen- 
tal value. 




more is that the magnetic moment comes out to be the 
experimental values at the same time. 

We have demonstrated that it is possible to perform 
highly precise calculation of the total energy in order to 
obtain both the correct easy axes and the magnitudes 
of MAE for Fe and Ni. This has been accomplished by 
including the strong correlation effects via taking intra- 
atomic repulsion and exchange into account, and incor- 
porating the non-coUinear magnetization. 

z 




y 



U(eV) 

FIG. 6: Brillouine Zone of fee crystal strueture. 



FIG. 5: The magneto-erystalline anisotropy energy MAE = 
£(111) - _B(001) functions of U. The experimental MAEs 
are marked by arrows for Fe (1.4 fxeV) and Ni (—2.8 fxeV). 
The values of exchange parameter J for every value of U are 
chosen to hold the magnetic moment of 0.61 /is in Ni and 2.2 
in Fe 

For Ni, we walked along the path of parameters U and 
J which hold the magnetic moment to 0.61 ^lb- The 
MAE first increases to 60 ^leV iJJ = 0.5 eV , J = 0.3 eV) 
and then decreases (see Fig. ||). While decreasing it 
makes a rather flat area from U = 1.4 eV , J = 0.9 eV 
to U = 1.7 eV, J = l.l eV where MAE is positive and 
around 10 iieV. After the flat area, the MAE changes 
from the wrong easy axis to the correct easy axis. The 
correct magnetic anisotropy is predicted at U — 1.9 eV 
and J — 1.2 eV. The change from the wrong easy 
axis to the correct easy axis occurs over the range of 
SU ^ 0.2eV, which is of the order of spin-orbit coupling 
constant O.le^). 

For Fe, the MAE is calculated along the path of U and 
J values where the magnetic moment is fixed to 2.2 /i^. 
At = and J = eV, the MAE is 0.5 fieV. It 
increases as we move along the contour in the direction of 
increasing U and J. The correct MAE with the correct 
direction of magnetic moment is predicted at ?7 = 1.2 eV 
and J = 0.8 eV. 

In this study of effects of U and J, it is again possi- 
ble to predict the correct magnetic anisotropy. What is 




FIG. 7: Calculated band structure of Ni at [/ = 1.9 eV and 
J = 1.2 eV. The solid and dotted lines correspond to majority 
and minority dominant spin carriers. 



IV. MAE AND BAND STRUCTURE 

We now present implications of our results on the cal- 
culated electronic structure for Ni. Fig. ^ shows Bril- 
louine zone and Fig. ^ shows band structure of Ni. We 
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L-neck 




FIG. 8: Calculated Fermi Surface of Ni at [/ = 1.9 eV and J = 
1.2 eV . The solid and dotted lines correspond to majority and 
minority dominant spin carriers. Dominant orbital characters 
are expressed. Both experimentally confirmed pocket and 
L neck can be seen. The X2 pocket is missing, which is in 
agreement with experiments. 



now discuss the band structure at the T point. From the 
below we see three majority bands and three minority 
bands. The lowest one is dominated by Asp orbital, the 
middle one is dominated by t2g orbital, and the highest 
one is dominated by Cg orbital. In this paragraph, we 
call the bands by its dominant orbital character at the T 
point. The eg band has two subbands. One of the two 
subbands is strongly hybridized with sp bands from the 
middle of the VX line. This band can be easily identified 
since the energy grows rapidly from the middle of the TX 
line and the rapidly growing part is dominated by sp or- 
bitals. Inspecting the sp band, we observe that the band 
is rather flat at later part of the VX line. In this area, 
the band is dominated by eg orbitals. The other subband 
of eg band remains unhybridizcd along the TX line. The 
t2g band has two subbands in the TX direction. One of 
the subbands is doubly degenerate and the other band is 
singlet. The doubly degenerate band can be identified by 
looking at the X point since this band breaks into two 
subbands there. Fig. || shows the Fermi surface deduced 
from the band structure. 

We now discuss the band structure in relation with 
magnetic anisotropy. One important feature which 
emerges from the calculation is the absence of the X2 
pocket (see Fig. ||). This has beenjiredicted by LDA but 
has not been found experimentallyEj. Inspecting the four 
band structures in band structure figures near the Fermi 
surface (Fig. ||), we notice that as we move away from 
the TX direction to TL direction, there is a band, which 
is above the Fermi surface at the X point, is submerged 
below the Fermi level. The band is the doubly degener- 
ate t2g subband and the pocket generated by this band 
is the experimentally confirmed X5 pocket. 

In the LDA band structure (see the top of Fig. |^), 
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FIG. 9: The X pockets. From the top, band structure near 
X of LDA, LDA+U &t U = 4, eV J = Q eV , LDA+U at 
U ^ l.<j eV J = 1.2 eV, and LDA+U at [/ = 1.7 eV J = 
1.1 eV . The size of the windows are 1 Ry. The solid lines 
represent bands with dominating up-spin. The dotted lines 
represent bands with dominating down-spin. X = (0,0,1), 
A = (1/16, 1/16, 1), B = (2/16, 2/16, 1), C = (3/16, 3/16, 1). 



there is another band forming another X pocket. This 
is the unhybridized eg subband and the pocket gener- 
ated by this band is the LDA X^ pocket. The LDA X2 
pocket has not been found experimentally. In LDA-I-U 
(see the second and the third of Fig. |^), the correspond- 
ing band is pushed down below the Fermi level and no X2 
pocket is present conforming to the experiments. This is 
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expected since correlation effects are more important for 
slower electrons and the velocity near the pocket is rather 
small. For the path with J = 0, we find that the removal 
of X2 pocket is around U — 2> eV far off from U = 4 eV 
where the correct magnetic anisotropy is predicted. No- 
tice that the corresponding band is way below the Fermi 
level at U — A eV (see the second of Fig. |^). For the 
path with the experimental magnetic moments, we also 
find that the removal of the X2 point is near the point 
U = 1.9 eV and J = 1.2 eV where the correct magnetic 
anisotropy is predicted. A very important point to notice 
is that the band that makes X2 pocket in LDA is just be- 
low the Fermi level (see the third of Fig. This brings 
a suspicion that the point U ~ 1.9 eV and J = 1.2 eV 
is where the X2 pocket just disappear. For this reason 
we study the X2 pockets at C/ = 1.7 eV and J = 1.1 eV 
(see the bottom of Fig. The corresponding band is 
just above the Fermi level forming a tiny pocket . This 
removal of X2 pocket near the point where the correct 
magnetic anisotropy is predicted, strengthens the con- 
nection between MAE and the absence of X2 pocket. 
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FIG. 10: Diagonal Components of Self energy matrix E;^^ for 
minority spin. Notice that Cg band (circle) is more pushed 
down than t2g band (square) 

This disappearance of X2 pocket can be explained in 
terms of self energy . For simplicity, we consider 
only diagonal components. Due to the cubic symmetry, 
these diagonal components are the same for orbitals in 
t2g band and for orbitals in eg band: 



^aa — 



t2g 



for orbitals in t2g bands 
for orbitals in eg bands 



(37) 



Since both X2 pocket and X^ bands are dominated by 
minority spins, we consider only minority spins here (see 
Fig. |l^. At U = 1.9 eV and J = 1.2 eV, we find that 
Et, = -0.09 eV and Ee = -0.30 eV. This shows that 
Cg bands are suppressed more by strong correlation than 
t2g bands. Since X2 pocket is dominated by eg bands 
and X5 pocket is dominated by t2g bands, the X2 pocket 
is removed by Strong correlation effects while X^ pocket 
survives. For comparison, we find that St^^ — —0.04 eV 
and Ee = -0.15 eV at U = 1.1 eV and J = 0.64 eV 

There has been some suspicions that the incorrect po- 
sition of the X2 band within LDA was responsible for the 
incorrect prediction of the easy axis within this theory. 



Daalderop and coworkerS removed the X2 pocket by 
increasing the number of valence electrons and found the 
correct easy axis. We therefore conclude that the absence 
of the pocket is one of the central elements in determining 
the magnetic anisotropy, and there is no need for any ad- 
hoc adjustment within a theory which takes into account 
the correlations. 

We now describe the effects originated from (near) de- 
generate states close to the Fermi surface. These, have 
been of primary interest in past analytic studieaS'Ll. We 
will call such states degenerate Fermi surface crossing 
(DFSC) states. The contribution to MAE by non-DFSC 
states comes from the fourth order perturbation. Hence 
it is of the order of A*, where A is spin-orbit coupling 
constant. The energy splitting between DFSC states due 
to spin-orbit coupling is of the order of A because the 
contribution comes from the first order perturbation. Us- 
ing linear approximation of the dispersion relation e(fcA), 
the relevant volume in fc-space was found of the order 
A'^. Thus, these DFSC states make contribution of the 
order of A'*. Moreover, there may be accidentally near 
DFSC states appearing along a line on the Fermi sur- 
face, rather than at a point. We have found this to be 
the case in our LDA calculation for Ni. Therefore the 
contribution of DFSC states is as important as the bulk 
non-DFSC states though the degeneracies occur only in 
small portion of the Brillouine zone. 

This importance of the DFSC states leads us to com- 
parative analysis of the LDA and LDA-f U band struc- 
tures near the Fermi level. In LDA (see top of Fig. I), 
five bands are crossing the Fermi level at nearly the same 
points along the FX direction. Two of the five bands 
are degenerate for the residual symmetry and the other 
three bands accidentally cross the Fermi surface at nearly 
the same points. There are two sp bands with spin up 
and spin down, respectively. The other three bands are 
dominated by d orbitals. In LDA+U (see second and 
third of Fig. ^), one of the d bands is pushed down be- 
low the Fermi surface. The other four bands are divided 
into two degenerate pieces at the Fermi level (see Fig. 
^: Two symmetry related degenerate di bands and two 
near degenerate sp^ and spi bands. In sum the corre- 
lation weakens the effect of degenerate bands along FX 
direction. 

In LDA (see the top of Fig. we found that two 
bands are accidentally near degenerate along the line on 
the Fermi surface within the plane TXL. One band is 
dominated by d| orbitals. The other is dominated by 
orbitals near X and by S| orbitals off X. This accidental 
DFSC states persist from FX direction to FL direction. 
(Along FL direction, the degeneracy is for the residual 
symmetry). In LDA-I-U, these accidental DFSC states 
disappear(see Fig. ||). With the correlation effect U, this 
accidental DECS states along a line on the Fermi sur- 
face move away from the Fermi surface leaving only the 
states along FL direction DFSC. This degenerate states' 
moving away from Fermi surface makes the first order 
perturbation effect sums up to zero. 
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FIG. 11: The L neck of LDA+U at [/ = 1.9 eV J = 
1.2 eV . Notice the absence of L pockets, nspecting the 
figures, we see that there is a L neck. The solid lines rep- 
resent bands with dominating up-spin. The dotted lines 
represent bands with dominating down-spin. The hori- 
zontal line is the Fermi level. The points are: A'^ = 
(29/64,29/64,19/32), O = (30/64,30/64,18/32), P = 
(31/64,31/64,17/32, Q = (125/256,125/256,67/128), 
(126/256,126/256,66/128), 5" = (127/256,127/256,65/128), 
L = (128/256, 128/256, 64/128). 



We now discuss another feature of Fermi surface, that 
is the L neck. The experimental L neck is a spin-up 
dominated band. From the Fig. [l|, we see one spm- 
up dominated band is just below the Fermi level in TL 
direction. As we move away from the TL direction, we 
see that the band is surfacing up above the Fermi level. 
This is the experimentally confirmed L neck. This L neck 
can be found at both (f/, J) = (1.9 eV, J = 1.2 eV) and 
(t/, J) = (4.0 ey, J = eV). For (C7, J) = (4.0 eV, J = 
eV), we find additional L pockets that has not been 
found in experiments. 

Based on the tight-binding.j|nodel, the importance of 
DFSC states has been showrolll. We see that strong cor- 
relations reduce number of DFSC states in TX direction 
and remove the near degenerate states on VXL plane. 
We conclude that the change of DFSC states is another 
important element that determines the easy axis of Ni. 
The correct configuration is tantamount to the correct 
Fermi surface. We see that it is possible to find the cor- 
rect Fermi surface at f/ = 1.9 eV and J = 1.2 eV where 
the magnetic moment and magnetic anisotropy can be 
correctly predicted. A point to note is that the Fermi 
surface does not come out right for other points in ([/, J) 
parameter space. 
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FIG. 12: The band structure of Fe for LDA. 
Now we discuss band structure of Fe. Unlike Ni, strong 
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FIG. 13: The band structure of Fe for LDA-I-U at t/ = 1.2 eV 
and J = 0.8 eV. 



correlation effects do not change bands structure of Fe 
near the Fermi surface significantly (see Fig. |l^ and ^3|). 
Most of the minority spin d bands lie above the Fermi 
surface and the self energy is not large enough to push 
the bands down near the Fermi surface. This is very dif- 
ferent situation from that of Ni, where minority spin d 
bands are near the Fermi surface and the self energy is of 
comparable magnitude. The self energy of minority spin 
eg band is positive for Fe. Since most of this band lies 
above the Fermi surface, the effects of strong correlation 
to occupied states is insignificant. Majority spin d bands 
lie far below the Fermi level. Since the self energy for 
this band is negative, the bands are pushed down more 
by strong correlations. This is the same situation as that 
of Ni and there is no significant changes near the Fermi 
level by strong correlations. The band structure of LDA 
and LDA+U are drawn in Fig. |l2| and |l^ . Compar- 
ing the two band structure, one can see that there is no 
significant change near the Fermi surface. This explains 
why magnetic anisotropy of Fe has been solved in LSDA 
while magnetic anisotropy of Ni needs a proper treatment 
of strong correlations. 



V. SPIN-OTHER-ORBIT COUPLING 

Our implementation of relativistic effects neglects 
other relativistic corrections than spin-orbit coupling. In 
this section, we first present that other relativistic effect 
can be seamlessly incorporated into out approach, then 
discuss spin-other-orbit coupling, which might be impor- 
tant in calculating magnetic anisotropy energy. 

We can divide relativistic corrections into one-body in- 
teractions and two-body interactions. One-body interac- 
tions are spin-orbit coupling, relativistic mass correction, 
and Darwin term. These one-body interactions can be in- 
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eluded into Kohn-Sham Hamiltonian in the same way as 
spin-orbit coupling is included. Only spin-orbit coupling 
terms is relevant to magnetic anisotropy energy and it 
has been the main theme of this paper. 

There are two two-body interactions coming from rel- 
ativistic corrections, spin-other-orbit coupling and spin- 
spin interaction: 



Hsoo - 



1 

^E^^"«J ^P»-(^'' + 2«^)' (38) 



2 ^ 



3, 



.(39) 



We now present the scheme of including two body inter- 
actions. A fermionic field operator 'J can be expanded 
along orthonormal complete bases {^a} where a is a 
combined index: 



(40) 



In this bases, the canonical anti-commutator relation can 
be written down in creation and annihilation operators: 
{ajj,a/3} — 5ap- In second quantization approach, any 
two particle interaction can be written down as 



Hint = ^ X! "'WjasapU 



(41) 



The matrix element Ua^s^yS can be found by 



af3j5 



$i(ri)$t(r2)„(l,2)$5(r2)$/3(ri)dridr2 



(44) 



(42) 

where v{l, 2) is a two particle interaction between parti- 
cle 1 and particle 2. The expectation value of Hi^t can 
expressed as 

(i?int) = - na,sn^f3)U"^^^ . (43) 

a/375 

Introducing another interaction matrix 

Ja/?75- J $i(ri)$;(r2)f(l,2)$^(r2)$5(ri)dridr: 

The Eq. ^ can be written as 

(Hint) = - J^'^'^'Wan^s- (45) 

The approach described in section ^ can be straight- 
forwardly extended to accommodate general two particle 
interactions by allowing spin dependent interaction ma- 
trices [/"'^T^ and J"/37<5. Therefore, all the relativistic 
corrections can be included in the framework of DMFT 
and its static limit LDA-I-U. 

The Coulomb interaction we have considered is one of 
the two particle interaction. In case of Coulomb interac- 
tion the interaction matrix Uap-yp is independent of spin. 



For more rigorous discussion we decompose the index a 
to (a,cr), where a is spin index and a contains all the 
others. The interaction matrix of Coulomb interaction is 



UaP^P — Uabcd^aa(7b^<^cf7d 



(46) 
(47) 



The Uabcd and Jabcd interaction matrices are the same 
interaction matrices in Eq. |l^. 

For general two particle interaction, the self energy is 
for spin-diagonal elements 





•^abcd I'^cd 


cd 




+ Y^Kbcf - 


' ■-'abed J^cd 


cd 






TCTCTCT-CTN cr-cr 
•^abcd /"-cd 


cd 




+ Y^Kb^r 


Taa-aa\ -aa 
■^abcd )"'cd 



(48) 



cd 

Unlike the case of Coulomb interaction, the spin- 
diagonal elements of self energy can get contributions 
from spin-off-diagonal density matrix components. No- 
tice that there are contributions from spin-diagonal ele- 
ments with opposite spin through exchange interaction 
(Jabcd'^'^n^d'")- The spin off-diagonal elements of self 



energy S^j^'^ takes the same form as the spin-diagonal 
elements except that the spin corresponding to the index 
b is —a instead of cr. 

We now study spin-other-orbit coupling. In this 
interaction, the spin degrees of freedom sees spatial 
anisotropy through orbital degrees of freedom. Spin- 
other-orbit coupling, therefore, may be important in 
magnetic anisotropy, though its contribution to MAE 
might be small. Since each term in spin-other-orbit cou- 
pling contains only one spin operator, either spins of a 
and f3 or spins of 7 and S must be the same in Eq. 
Using this observation, the spin-diagonal elements of self 
energy is simplified to 



Kb - YiK^Td'^^'^^n-r" (49) 

cd 

I \^ [TT^^^^ j(Taaa\ aa 
2_^\^ahcd ~ 'Jahcd )^cd 

cd 

I \^ /jraaa-a _ Taaa-a\ a-cr 
Z_^\^abcd '^abcd )"'cd 

cd 

2_^\^abcd ^ '^abcd )^cd 
cd 

Notice that the contribution from spin-diagonal density 
matrix components has the same form as the one of 
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Coulomb interaction except that U"^^J^ need not to be 
the same as Uabcd'^~'^ ■ Since the spin-off-diagonal den- 
sity matrix elements are small compared to the spin- 
diagonal elements, we may neglect such contributions. 
When Coulomb interaction is consider as well as spin- 
other interaction, its contribution to the interaction ma- 
trix elements of the form f/^r , Kbd'^"' ^ J--ff is 
much larger than that of spin-other-coupling. In this 
approximation, the contribution of spin-other-orbit cou- 
pling is negligible and the spin-diagonal elements of self 
energy is then the same as that of Coulomb interaction. 
Spin off-diagonal elements is given by 

^ab ~ Z_^\ '^abcd )"'cd 
cd 

+ 2^(C^a6cd -J abed hcd (50) 

cd 

I \ ^ ( J ja — aaa _ tct — crcrcr\ aa 
^ 2_^\^abcd ~ '\bcd J^cd 
cd 

The contributions from spin-off-diagonal density ma- 
trix elements has the same from as that of Coulomb 
interaction. Notice that there are contributions from 
spin-diagonal density matrix elements. This contribu- 
tion should be treated on the same footing as the contri- 
butions from spin-off-diagonal density matrix elements, 
hence the correction from spin-other-orbit coupling is 
not negligible. For an extensive study of this effect, we 
need an implementable formula. We developed an for- 
mula that can be implemented in Appendix 

An extensive study on the effects of spin-other- 
coupling to spin-off-diagonal elements of self energy re- 
quires another paper. Instead, we estimate the order of 
magnitude of spin-other-orbit coupling here. The order 
of magnitude of spin-other-coupling can be estimated by 
using (l/rfj) - (1/rijf to {C I / 2amc^f C I , where CI is 
the strength of Coulomb interaction and a is the fine 
structure constant. We see that the order of magnitude 
of CI is 1 eV for d shell electrons of Ni and Fe. Therefore, 
the spin-other-interaction is 10~® times smaller than the 
Coulomb interaction. The effect of spin-other-orbit cou- 
pling is, therefore negligible. This-tas been confirmed by 
the study of Stiles and coworkerscfl 



VI. CONCLUSION 

We have studied the effects of strong correlations on 
magnetic anisotropy. We find that magnetic anisotropy 
changes as correlation strength U and J change. These 
has been shown along several paths in ([/, J) space. The 
correct magnetic anisotropy can be predicted at several 
points in ([/, J) space for both Ni and Fe. Moreover, 
magnetic moment, magnetic anisotropy, and Fermi sur- 
face are correctly predicted simultaneously at physically 



acceptable values of U and J: U = 1.9 eV and J = 1.2 eV 
for Ni and U = 1.2 eV and J = 0.8 eV for Fe. 

It is remarkable that the values of U necessary to re- 
produce the correct magnetic anisotropy energy are very 
close (within 1.2 eV) to the values which are neededjto de- 
scribe photoemission spectra of these materiala^j'LJ. The 
correct estimation of U is indeed a serious problem. Dif- 
ferent estimates based on different methods give a range 
of values from 1 eV to 6 eV. The work of Katsnelson 
and coworkers used a dynamic approach utilizing finite 
temperature Quantum Monte Carlo method whereas we 
use a static approach at zero temperature. Considering 
these discrepancies, an agreement within 1.2eV shows an 
internal consistency of our approach and emphasizes the 
importance of correlations. 

For Ni, we find that the LDA X2 pocket disappears 
near U = 1.9 eV and J = 1.2 eV where the correct 
magnetic moment and the correct magnetic anisotropy 
are predicted simultaneously. This suggests a tight re- 
lation between the absence of X2 pocket and the cor- 
rect magnetic anisotropy as suggested by previous work. 
The change of DFSC states from LDA to LDA-I-U is dis- 
cussed. DFSC states are old suspects of dominant con- 
tribution to magnetic anisotropy. The correct configura- 
tion of DFSC states and the absence of X2 pocket are 
by-products of the correct Fermi surface that we predict 
for the first time. For Fe, we find that strong correlations 
do not change band structures significantly. This is at- 
tributed to the fact that most of minority bands lie above 
the Fermi level, hence the change of band structures has 
no large effect to occupied states. This distinguishes the 
magnetic anisotropy problem for Fe from the one for Ni. 

The calculations performed are state of the art in 
what can currently be achieved for realistic treatments 
of correlated solids. Despite the great successes of the 
LDA-fU theory in predicting material properties-ef corre- 
lated solids (for a review, see book of AnisimovEJ), there 
are obvious problems of this approach when applied to 
metals or to systems where the orbital symmetries are 
not broken. The most noticeable is that it only describes 
spectra which has Hubbard bands. A correct treatment 
of the electronic structure of strongly correlated electron 
systems has to treat both Hubbard bands and quasipar- 
ticle bands on the same footing. Another problem occurs 
in the paramagnetic phase of Mott insulators, in the ab- 
sence of any broken symmetry the LDA -I- U method 
reduces to the LDA, and the gap collapses. In systems 
like NiO where the gap is of the order of eV, but the 
Neel temperature is a few hundred Kelvin, it is unphys- 
ical to assume that the gap and the magnetic ordering 
are related. For this reason the LDA-t-U predicts mag- 
netic ordepin cases that it is not observed, as, e.g., in the 
case of PrO. Since LDA-f U method is a limiting case of 
more general DMFT, these difficulties are expected to be 
overcome in DMFT. Further studies should be devoted 
to improving the quality of the solution of the impurity 
model within DMFT and extending the calculation to 
finite temperatures. 
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APPENDIX 

In this appendix, we develop an computationally im- 
plementable formula for spin-other-orbit coupling: 



v{l, 2) = ifj^B 3 

''12 



(tJi 2ct2) + 1 ^ 2, (51) 



where ai applies to states at ri and a2 to states at r2. 
For simplicity, we will only work the first term out. The 
formula for the second term can be found by exchanging 
1 and 2. Using the identity 



^ d r X L 
V = r— - I — 5—, 
or r'^ 

the spin-other-interaction can be written down as 

W(l,2) = -^I^Ll- ((71-^2(72) 
'^12 

n (fl) X f 2 . „ X 9 



(52) 



(53) 



dri 



1 12 

2 ^2 f 2 • f 1 - , . 
- i^B- • (tTl -f 2(72) 

'1 '12 

+ Ms 3 • ((71 -I- 2(72)- 

We will write the interaction matrix as 

= C/„^^5 -I- f7^^^5 -I- (54) 

where U^^^g comes from the first term, U'^^^g from the 
second term, and U^^^g from the third and fourth terms 
of the Eq. ^ For easiness of following index structure, 
we define the interaction matrix U as 

Uaisjs^ J $Ui-i)*]3(r2)w(l,2)$^(r2)$5(ri)dridr2 



We will consider bases in the from of 



(55) 



(56) 



where Xcr = (1, 0)"^, or (0, 1)"^. Note that the angles d and 
(j) is measured with respect to the origin of a coordinate 
system. 



We want to separate the spin-other-orbit interaction 
into radial part and angular part. This can be done by 
using the following property 

^ =E/'(^l'''2)>'i™(^2,02)r;„(01,(/.l), (57) 



where 



fi{ri,r2) = 47r 



l_i {rf + — 2rir2.T)2/3 
For example we write down /o and /i explicitly 

fo[ri,r2) = [-. r ■ 

rir2 \ri - r2\ n + r2 



dx. (58) 



(59) 



fiiri,r2) 



2tt (ri - r2f +rl + rl 



Vi - 7-2 1 
{ri + r2f +rl+rl ^ 
n + r2 



(60) 



With these identities, we can calculate a readily imple- 
mentable formula for spin-other-orbit interaction. After 
some manipulation, we can write the contribution from 
the first terms of the Eq. M as 



uh 



{iilimicn){i2hm2cr2)(i3hm3a3){i4,Um4a4) 

X 



= -A^lX! / dridr2fi{ri,r2) 

Im 



{l2m2\Yim\l3m3){limi\YimOcria2a3aAkm4:), (61) 

where the operator Ol.^„^^^^^ is defined by 

2S...Aa3+S.2.30l,^, (62) 



^aiCT2(T3Cr4 



— Lz 

oil = 

0]i = 

= L+ 

The contribution from the second terms of the Eq. |5i 
reads 



U, 



(il/imiCTi)(i2'2m20'2)(i3'3™3 0-3)(i4'4"l4 0-4) 

fl,f2) X 



Im 



d 



'-'^/imicri)(i2ni20-2)(/3m3<T3)(i4m40-4)(^'^)' (^^) 



where 



^(/imicri)(;2m2 0'2)(i3™30'3)('4m40'4) (^'^) 

{Im) 

(Tl(T4 



CT2CT3 (iimi)(i2m2)((3m3)('4m4) ^ ^ ' 
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and the operator N??-'^'' w, Aim) is defined 



27r 



as 



N, 



TT 

(iimi)(i2m2)(i3m3)('4m4) 



(Im) 
(Im) 



(Iimi)(l27yi2)(l3m2,)(limi) 
^ (Iimi)(l2m2)(l3m3)(limi) i^^) 
= ^^thmi){l2m2){hm3)(Um4)i^'^) 

(Im) 

^ ^(iimi)(;2m2)(i3m3)(i4m4) 

(Im) 



{Iimi){l2m2){hm3)(l4,m4) 



{Iimi){l2m2){l3m3)(l4,m4) 



A 



{hmi){l2m2){l3m3){l4m4) 



(Im) 



(65) 
(66) 
(67) 
(68) 



The operator A(^i^rni){i2m2)ihm3){i4m4){lm) is defined by 
a cross product of another operator 5((imi)(;2ni2)- 

^ihnii){l2m2){l3m3)iUni4)ii'm) 

= ■S(;imO(i4m4)(M X ■B(i2m2)(i3™3)(^^")' (69) 

where 



BLi'm'ilrn) = ^/y(?m|r;,„(ri_i -yn)|;'m'X70) 
2^ 



^Lr™'(M = ^/y(/m|rz™(ri_i +yn)|/'m'X7i) 
BLi'm'ilm) = — {lm\YirnY,o\l'm') (72) 
The contribution from the last two terms of the Eq. M 



-O 



(iimiCri)(i2m2<T2)('3™3f3)(Um40'4 



i(/m) 



(73) 



The operator O?, ^ „ \n ^ \n ^ ^ mi ^ ^ ^('"1)7 which 

^ {Limiai){l2m2<72)\h^3f^3}[^4m4a4)^ /' 

represents the contribution from the third term, is de- 
fined by 

n3 (Im) 

^(limicri)(l2m2<J2){l3m3<73){l4m4a4) V / 
= {l2m2\YimYiQ\l3m3) X 

{hmi\YimYioOl^„^„^^^ \Um4) 

-{l2m2\YlmYii\l3m3) X 



(Zimi \YimYi-iOl^„^„^„JUm4) 

-{l2m2\YijnYi^l\l3m3) X 

{hmilYimY^iOl^ 



where 0\,^„„„„, is given by the Eq. 62. The opera- 



''■2<''3''"4 



1/47714 



ClO'20"3^7"4 

tor Of, X,, X,, s,, s (/m) which represents 

(iimi(Tl)(i2m2<T2)(i3m3<T3)(i4rri4CT4) ' ^ 

the contribution from the last term, is defined by 

dm) 

^{Iimiai)(l2m2cr2){l3m3a3){l4m4a4) V'"-) 

= y/2{l2m2\Yi^Yio\hm3) x 
{limi\YimOl^„^„^„^L^\hmi) 

-(^2m2|yimYll|^3"^3) X 

(/imi |yi™0^^^^^3^^L_ 1/4^4) 
-{h■m2\Ylr,^Yl^l\l3m■i) x 
{hm,\YunOl^,^_,^,^L+\U7n^), (75) 

where 0^i(T20-3<j4 is given by the following equations: 



^fTiCr2Cr3a-4 ^'^^7'lCJ'4^(J2a-3 



'5^2.3 0^, .4 (76) 



(ii;inii(Ti)(i2i2m2 0"2)(i3'3m30'3)(i4'4"l40'4) 
<;i(^l)'^»2i2(''l)'?^«3/3('^l)'^»4/4(''l) X 



47r 



-O; 



(iimi(Ti)(i2ni2 0-2)('3m3'T3)('4"l40"4) 



(Im) 





= \/2Yio 




= -V2rio 








= -1^11. 
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